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Abstract 

Stellar astrospheres — the plasma cocoons carved out of the interstel- 
lar medium by stellar winds — are continually influenced by their passage 
through the fluctuating interstellar medium (ISM). Inside dense interstellar 
regions, an astrosphere may be compressed to a size smaller than the liquid- 
water habitable zone distance. Habitable planets then enjoy no astrospheric 
buffering from the full flux of Galactic cosmic rays and interstellar dust and 
gas, a situation we call "descreening." Recent papers (Yeghikyan and Fahr, 
Pavlov et al.) have suggested such global consequences as severe ozone de- 
pletion and glaciation. Using a ram-pressure balance model that includes 
gravitational focusing of the interstellar flow, we compute the size of the as- 
trosphere in the apex direction as a function of parent star mass. We derive 
a dependence on the parent-star mass M due to gravitational focusing for 
densities larger than about 100 {M/Mq)^^ cm^^. We calculate the inter- 
stellar densities required to descreen planets in the habitable zone of solar- 
and subsolar-mass stars and find a critical descreening density of roughly 
600 (M/Mq)"2 cm"^ for the Sun's velocity relative to the local ISM. Fi- 
nally, we estimate from ISM observations the frequency of descreening en- 
counters as 1-10 Gyr^^ for solar-type stars and 10^-10^ times smaller for M 
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stars. Given this disparity, we conclude that M star habitable-zone planets 
are virtually never exposed to the severe effects discussed by Yeghikyan and 
Fahr and Pavlov et al. 



1 Introduction 



The heliosphere is the magnetic bubble carved out of the interstellar medium by 
the expanding solar wind. A recurring idea is that during the Sun's lifetime en- 
counters with dense interstellar environments have compressed the heliosphere, 
enhancing cosmic ray fluxes and hydrogen accretion, affecting atmospheric chem- 
istry (inc luding o zone), and perhaps altering s urface mutation rates and c 
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20081) . There may even be em- 



pirical evide nce for such an encounter in lunar soil samples (IWimmer-Schweingruber and Bochslei 
2000l l200lb . A similar variation should appear for other stars, given that "astro- 
spheres" around other stars have no w been directly observed through their inter- 
action with the interstellar medium (" Wood et a/.Ll2002Ll2005b . 

Interest in the potential habitability of planets orbiting stars of mass smaller 
than the Sun's — cool, red, faint stars of spectral class M with masses between 
0.08 and 0.6 Mq — has increased in recent years. For simplicity, we refer to all 
subsolar-mass stars as M stars, since K main-sequence stars cover such a small 
range in mass. The liquid- water habitable zone (H Z) is defined e ssentially by a 
certain range in received flux at a planet's surface ( Kasting et all IT993i) . and M 
star planets have small luminosities, so the habitable-zone distance rnz can be 
extremely small, from 0.3 AU to 0.02 AU as mass varies from 0.6 to 0.08 Mq. 
Habitable planets around M stars may be common for many reasons, including 
the large number of M stars compared to any other class of star, the detection of 
extrasolar planets orbiting M stars, their presumed ability to form terrestrial-mass 
planets if they have sufficiently massive disks, their obser ved disk frequency , and 
even evideiice for d ebris disk substructure in one case. See lScalo et a/.l(|2007b and 
Tarter et all (|2007b for detailed reviews of these concepts. 

The character of the astrospheric environment expected for habitable-zone 
planets orbiting stars of different masses has not been examined. The present work 
focuses on the frequency of reductions in the size of a star's astrosphere to within 
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the HZ or smaller induced by passage through the fluctuating conditions in our 
Galaxy's interstellar medium (ISM). Such astrosphere size changes are inevitable 
given the range of conditions in the ISM and may result in significant changes in 
the Galactic cosmic ray (GCR) flux, as scattering of cosmic rays in the astrosphere 
is the main filter (besides any possible planetary magneto sphere) shielding plan- 
ets from G CRs. The effects of GCRs on planetary atmosph eres are significant and 



varied (e.g. iGriessmeier et all |2005|; iGrenfell et all |2006|) . In addition, changes 



in astrosphere size should lead to large variations in the planetary accretion rate 
of interstellar hydrogen or dus t, a possibilit y that has been discusse d for several 



decades (Fahr, 1968; MccSjl975; Bzow ski et all [l996l: IVa bush ita and AUenl 



19891 : fplorinski et alii 120031: lYeghikvan and Fahil l2004bl: IPavlov etall l2005bh : 



the climatic changes associated with such accretion could be varied and severe 

I II — — I I — — II — 

(e.g. the noctilucent cl oud scenario of [M cKay and Tho maslll978 ': McKayl ll985 . 

or ozone depletion as inl Yeghikyan and Fahr.,2004b. : Pavlov et a/.li2005b,) . 

In the present paper we show that the size of astrospheres around M stars will 

be controlled by the star's gravity during descreening episodes and estimate the 

frequency of these episodes for M stars as compared to solar-type stars. Section |2] 

will describe the analytic models for astrospheres that we will use. Section |3]gives 

the resulting mass and age dependences, and Section [33] will use these models to 

estimate the frequency of descreening encounters for solar- and subsolar-mass 

stars. 



2 Methods 

2.1 Astrosphere Model 



Parked d 1963b showed that the location of the termination shock of the solar wind 



given by the shock-jump conditions was roughly where the ram pressure of the 
outflowing solar wind balanced the pressure of the counterstreaming ISM: 

Asm = ^ ^ , (1) 
r 

where is the termination shock distance, ro is a reference distance, po is the 
solar wind mass density at ro, vq is the solar wind velocity at ro, and Pism is 
the total pressure of the ISM, including ram, thermal, and magnetic pressures. 
When the ISM is streaming supersonically and superalfvenically past the stellar 
system (such as for the present day solar system), its ram pressure dominates 
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i^SM, such that P^sm — Pism^^, where PisM is the total ISM density and V is the 
relative velocity between the star and the local ISM flow. Taking the size of the 
astro sphere to be the distance from the star to the termination shock, we have 



ra _ ^noy/^vo 
ro~\NJ V' 



(2) 



where n and A'^ will henceforth refer to the total ion number density of the solar 
wind and ISM, respectively. 

To estimate the stellar wind pressure we assume that the wind within astro- 
spheres of stars of different spectral types (or equivalently masses) are analogues 
of the solar wind. Little information on the wind is available for stars besides the 
Sun other than the variation of mass-loss rate with age and type for a few nearby 
stars dWood et a/.Ll2002ll2005|) . The wind speed vq in our model is taken to be 400 
km s~^ and assumed to be independent of di stance from the star, as is roughly the 
case for the solar wind in the ecliptic plane dWhang et all l2003b at radii beyond 
a few times the critical point (~ 0.1 AU). Recent work suggests that the steady- 
state wdiidjnayjn fact be slower and even su bsonic and subalfvenic near to the 
star (lErkaev et all 120051: IPreusse et all |2005|) . in which case the interstellar ram 
press ure would be balanced by wind magnetic and thermal pressure near M star 
HZs. iKhodachenko et all (|2006|) . however, shows that the steadily higher coronal 
mass ejection rate on M stars may significantly affect the wind speed in the HZ. 
To examine the possibility of a slow wind near the HZs of the lowest mass stars, 
we include a model in our results in which the wind speed is taken to be 50 km 
s^^ instead of 400 km s~^. 

For simplicity, we lump the wind density and velocity dependence on stellar 
mass and age into a scaling function: 



4.5 Gyr 



/3 



(3) 



where a and /3 are allowed to vary. Thus the ram pressure balance equation for 
the astrosphere size can be written: 



ra 
ro 



1/2 

V' 



(4) 



The in dices a and jS are extreme ly uncertain as present, but studies of stellar mass 
loss by IWood et all (|2002l 120051) suggest that younger and lower- mass stars may 
have stronger winds, implying a < and /3 < 0. 
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Observations of IWood et al\ (|2005b suggest j8 — — 2.33, but^ using the same 
relation between mass -loss rate and X-ray flux a s I Wood et all dlOOSb. th e re- 
cent surveys of X-ray fluxes of dwarf M stars by Penz and Micelal ( 2008 ) sug- 
gest /3 = —1.8. The mass dependence of the mass-loss rate is more complicated, 
how ever, and ob servations of some M stars show a lower mass-loss rate than the 
Sun ( Wood et al. l|2005t). If this is accurate, then it may be possible for M stars to 
have a lower mass-loss rate than the Sun, in which case, they would be easier to 
descreen. We include this possibility in our results in an example in which a = 1. 

Though crude, a pure ram pressure balance approximation is surprisingly ac- 
curate at predicting the solar wind termination shock distance and generally yields 
heliosphere sizes in the apex direction wi thin 10-20% of detailed multifluid , mul- 
tidimensional cornputational mode l s (e.g.lZank and Frischl.ll9 99: Fah r et a/.L 12000 ; 
Wood et ali 120021: Iplorinski et a/1 l2004l : IZank et al.L 120061: [Muller eTall l2008h . 



The obvious weakness of this analytic model is that, being one-dimensional, it 
cannot reproduce the full three-dimensional structure of the global heliosphere. It 
can, however, predict the minimum ISM densities above which planets in the HZ 
can be descreened for at least part of their orbit. In one sense, this model may even 
overestimate the astrosphere size and hence the descreening frequency because it 
neglects charge exchange with interstellar neutrals. 

We have up to now ignored the thermal and magnetic pressure contributions 
in the ISM. Thermal pressure would be important (for the solar velocity with 
respect to the local standard of rest) for ISM densities below about 0.1 cm^, 
where temperatures are high, but we are interested in the high densities that could 
lead to descreening, and these regions are typically cold {T ~ 10-20 K, Ferrierd . 
liool) compared to stellar velocities. For thermal pressure to dominate to ram 
pressure, the following inequality must be satisfied: 



T > 3000 



V 



lOkms 



K, 



(5) 



where V is the star- ISM relative velocity. For the Sun, V — 26 km s^^, so it would 
need to enter a region with T > 20,000 K for the thermal pressure to compare to 
the ram pressure in determining the heliosphere size, and these regions usually 
have densities well below the regimes we are considering he re. For instance , the 
warm ISM has a density below ~ 0.5 cm^-^ and T > 6000 K (iFerrierd . 1200 ih . 

Magnetic pressure is likely to be more important, although without a detailed 
MHD simulation of the entire ISM structure spanning orders of magnitude in 
scale, we can make only a rough estimate of the effect. Zeeman observations 
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([Crutchefl, 1 19991 : iBourke et all 1200 Ir) can be expressed in the form of an upper 
limit to a 5-density relation: 



B{N) 



N 



10 cm 



3 

1/2 

^1 i"G 



< 10 cm-3 
> 10 cm-3 



(6) 



At densities below 10 cm^^ this simply gives a constant pressure B^/Sn equal 
to 4 X 10^^^ dyne cm^^. Since the ram pressure is 8 x 10^^^ ^^20 dyne cm^^, 
the magnetic field isn't competitive with ra m pressure until dens ities below 0.05 



cm ^ (e.g. in the Local Interstellar Cloud; Florinski et a/.ll2004j have suggested 



that the magnetic field is anomalously large here). At higher densities, the square- 
root scaling gives a pressure that is proportional to the density, similar to the ram 
pressure. In this case the coefficient of the magnetic pressure acts to effectively 
increase the ISM-star velocity in t he ram pressure , but the effect is small, even 
after multiplying by a small factor (' Crutchei] . [l99 9') to statistically account for the 
fact that only the line-of- sight component of B is detected by the Zeeman effect. 
This is because typical star-ISM velocities exceed the root-mean-square turbu- 
lent velocity in the ISM, which is in rough equipartition with the magnetic field. 
Given these considerations, we neglect the magnetic field for simplicity, with the 
understanding that we somewhat overestimate the size of the heliosphere at low 
densities, and that in more realistic ISM models the magnetic field will undergo 
large fluctuations, often uncorrected with density, that may lead to episodes of 
magnetically controlled astrosphere sizes. 



2.2 Gravitational Focusing of the Interstellar Flow 

As planetary systems traverse dense clouds, gravitational foc using will de termine 
inflowing particle trajectories. This was first pointed out by iFahrl (|1968r) for the 
case of coUisionless orbits of proto ns in the inner solar system and then later 
first applied in the present context by iBegelman and ReesI (|1976h . Previous work 
has either neglected gravitational focusing or included its effect on the star-ISM 
veloc ity but only in the fluid limit (e.g. .Talbot et all 1 19761 : iTalbot and Newmanl 
1977b . 

Here we assume the ISM flow is coUisionless. The mean free path for Coulomb 
collisions is ~ 10^^ A^^^ cm, where A^ is the local ISM number density, and this 
is larger than the astrosphere size in most cases of interest. This assumption 
is marginal at low densities, but at the high densities required for descreening, 
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where gravitational focusing dominates the accretion, we show below that the as- 
trosphere size decreases with ISM density faster than does the coUisional mean 
free path, so the coUisionless assumption is valid. 

Neglecting collisions simply means that we neglect fluid effects, such as gas 
pressure, and assume the ISM particles follow ballistic trajectories until they "col- 
lide" with the wind-ISM interface, exerting a pressure. The physics of this interac- 
tion is beyond the scope of this discussion, but ram-pressure balance approximates 
the standoff distance between the wind and ISM quite well, and this is our primary 
concern. 

A com plete solution to s tellar accretion in the coUisionless limit was derived 
by loanby and Camm (1957). The infall density can be expressed as an integral 



that must be evaluated numerically in the general case, but for our assumption of 
small gas thermal velocity dispersion, Danby and Camm find that the density in 
the apex direction is 

where Nq is the unperturbed ISM number density far from the star, p = V/o,q^ = 
2GM/o^r, and o is the thermal velocity dispersion of the gas. The ram pressure 
on the astrosphere in the apex direction will then be this density multiplied by 

o 2GM 

V'^ = V^ + , (8) 

r 

where the second term on the right accounts for gravitational acceleration. Bal- 
ancing the outward wind ram pressure with the inward ISM ram pressure as before 
yields a quartic polynomial for the astrosphere size (see §3.21) . Most of the mass 
dependence of our result comes from the importance of the gravitational focusing 
terms involving in Eq. |7] 



3 Results 

3.1 Descreening in the High-density Limit 

The minimum cloud density Nf for focusing to dominate can be obtained from 
the condition b > Va, where b = 2GM /V^ is the distance at which the star's grav- 
ity becomes important, and is the pressure-balance astrosphere size neglecting 
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focusing. This gives 



Nf=\40^iM,t)Vio 



M 



-2 



cm 



-3 



(9) 



where V20 is V in units of 20 kni s^. Thus in dense molecular clouds (A^ ~ 10^- 
10-^ or larger; e.g. iLarsonl Il98ll : IScalol [1985) the accretion flow will be mostly 
determined by gravitational focusing except possibly for the very lowest-mass 
stars. For refere nce, the Sun 's present speed relative to the local ISM flow is 
about 26 km s~i (IWitteLl2004h . 

Assuming that gravitational focusing dominates, the density at the apex of the 
astrosphere is 

where A^o is the ISM density unperturbed by the star's gravity. The incoming ve- 
locity (much larger than the star-gas velocity V in this limit) is V' ~ (2GM/r)^/^. 
Expressed in terms of the gravitational radius b, we get 



and 



V'-V 



1/2 



(11) 



(12) 



where b/r 1 in the gravitational focusing limit. Equating NV'^ to the stellar 
wind ram pressure then allows a simple solution for the astrosphere size of a plan- 
etary system when the flow outside the astrosphere is dominated by gravitational 
focusing. The result is 



— ^4 



or, substituting for b. 



[a 



y2 4 3 

2G3m3 



ro 



^{M,t)no 



-2 



(13) 



(14) 



Notice that in this gravitational focusing limit the astrosphere size increases with 
increasing stellar velocity V. This unexpected behavior occurs because the slower 
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the star is moving, the stronger is the effect of focusing, so that the ram pressure in 
this case scales like V^^. And since the astrosphere size varies as Nq^ in this limit, 
while the mean free path of the incoming gas to particle-particle collisions still 
scales linearly with A'q"^ the astrosphere size is decreasing faster (due to focusing) 
than the mean free path is, and the validity of the coUisionless approximation, 
as measured by the ratio of mean free path to heliosphere size, increases with 
increasing density. Finally, the inverse mass dependence of the astrosphere size 
means that the coUisionless limit becomes rapidly more accurate for lower mass 
stars. 



By the definition of the HZ (neglecting albedo variations), rnz = (^/^g 



1 1/2 



AU where Lq is the bolometric luminosity of the Sun. Although no single power 
law can accurately represent the mass-luminosity relation for low-mass stars (see 
below), an adequate relation for our purposes is L/Lq = (M/M©)", with a ~ 
2-3 f or masses bet ween 0.1 and 1 Mp), based o n masses from suitable binary sys- 
tems ( Henry et al. I [1999; Del fosse et ali |2000() . Hipparcos distances, and empiri- 



cal bolometric corrections. The habitable-zone distance is then rnz = (M/M©)"/^ 
AU. Taking ro = 1 AU and the average number density of the solar wind at 1 AU, 
HQ, to be 7 cm^-^ and evaluating numerical factors: 

- 3.6xlo'»^)]!^^v"°'^ (15) 



The minimum ISM density for this descreening condition to hold is 

Nd>eOOV2o^{M,t)\y—\ cm-3, (16) 

where the numerical coefficient is approximate because the star-gas relative ve- 
locity will vary from star to star by factors of a few. For a = 2.5, the exponent of 
M is —2. 1, so a star with M = 0.5 Mq (early M star) will require an ISM density 4 
times larger than the Earth for complete descreening to occur, while a 0. 1 Mq star 
will require a density about 120 times larger. Our charge exchange model would 
yield an astrosphere a factor of (fo/V)^/^ smaller for a given ISM density, so it 
would yield a descreening density that is {vq/vY^^ times smaller than what the 
ram pressure approach above would yield. For any reasonable wind speed and 
ISM-star relative velocity, this would be only a factor of 2-3 at most. 

This result is easily understood as a combination of several effects. First, the 
smaller stellar mass produces less gravitational focusing, reducing the exterior 
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ram pressure and increasing the astrosphere size compared to larger masses. Sec- 
ond, the habitable-zone distance decreases with decreasing stellar mass, forcing 
the astrosphere to shrink to smaller distances for descreening of habitable plan- 
ets to occur. Third, the mass- and age-dependent wind scaling factor ^ (Mj) at a 
given age probably will increase with decreasing mass, leading to a stronger wind 
and larger astrosphere. 

So far we have assumed the limiting case in which the density field is com- 
pletely dominated by gravitational focusing in order to get an analytic result. In 
the following section, we solve the full quartic polynomial equation for the den- 
sity and give the exact (in the limit of a cold ISM) solution for ra/rnz and the 
critical density for descreening. 



3.2 Descreening in the General Case 

Following the notation of the previous section, the full polynomial equation for the 
astrosphere size can be derived in the coUisionless limit by equating the wind ram 
pressure with the ISM ram pressure, using the Danby and Camm (1957) gravi- 
tational focusing solution for the density at the astrospheric apex point and the 
velocity as modified by gravitational focusing (Eq. [8]). The result is 



1 



NV^ 



0, 



(17) 



where b is the gravitational radius as defined earlier and Pw{M, t) = t, (M, t) riQV^r^. 
To solve this equation, first make the substitution ra = u — b/2. This results in the 
depressed quartic equation: 



4 2 

u u — 



NV^ 



0. 



(18) 



Solving for u and taking the positive, real root as the astrosphere size: 




NV2 



l/2\ 1/2 



b 
2' 



(19) 



The quantity r^/rHz as a function of local ISM density is shown for three repre- 
sentative stellar masses (0.1, 0.5, and 1 Mq) in Fig. [T] 
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Figure 1 : Astrosphere size, in units of liquid- water habitable-zone distance (de- 
fined here as the distance receiving the same bolometric stellar flux as Earth), as a 
function of ambient interstellar gas density for stars of three different masses. The 
regime in which the inflowing interstellar ram pressure is dominated by gravita- 
tional focusing is at densities higher than the "knee" in each curve. The star- ISM 
relative velocity was taken to be 10 km s^^. The stellar wind speed was taken to 
be 400 km^^ in the first three cases and 50 km s^^ in the last case ^ representing 



the possibility of subsonic , subalfvenic winds in the HZ of M stars (lErkaev et at 



20051 : iPreusse et all 120051) . A slower wind speed seems to have a significant ef- 



fect on the astrosphere size, but even in this extreme case the astrosphere size is 
larger for the M star than for a typical solar-mass star for all ISM densities. 
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Figure \T\ shows that the astrosphere is much farther away in terms of the 
habitable-zone distance rnz for systems with lower-mass parent stars. The tran- 
sition into the gravitationally dominated regime is located, for each of the four 
model planetary systems, by the knee in the curves as the ISM density increases. 
This knee occurs for lower ISM densities for systems with more massive parent 
stars, since the effect of the star's gravity on the pressure balance grows at a given 
ISM density as mass increases. One case is shown with a slow, 50 km^^ wind 
velocity. This case shows that if M stars do have weaker winds than solar-type 
stars, then they will have correspondingly smaller astrosphere sizes. But even in 
the severe case of a wind speed of 50 km shown in Fig. [H the astrosphere is 
still larger relative to the HZ distance than for the solar-type stars. 

Figure [T] also shows clearly that an M star must encounter a very high density 
interstellar region in order to compress its astrosphere significantly, and therefore 
full exposures to the Galactic cosmic ray and accretion flux will be rare com- 
pared to a solar-type star. This occurs because denser interstellar regions occupy 
a smaller volume fraction of the Galaxy than lower-density regions, as discussed 
inCT 

A more quantitative calculation of the effect of stellar mass on the critical 
descreening density is obtained by solving for the density at which the astrosphere 
apex distance shrinks to the habitable-zone distance rnz, and then relating rnz to 
the stellar mass. From Eq. [l7]the critical descreening density is 



Nd = '-^^^f^[r-BZ + ^bri^+-b'rt,^ + -b\nzj ■ (20) 

The habitable-zone distance depends on the luminosity of the parent star, rnz ~ 
L^/^ in Earth-Sun units (again ignoring variations in albedo). Because the pri- 
mary (approximately constant) quantity that governs the evolution of the star is its 
mass, it is traditional to parameterize L in terms of stellar mass M. We used the 
data for stellar masses a nd lum inosities for stars less massive than 1 Mq given in 



-1/2 



Hillenbrand and Whitd (|2004|) and fit the following polynomial to it: 



logL = 4. lOlog^ M + 8. 16 log^M + 7.11 logM + 0.065, (21) 

where L and M are in solar units. The slope of this log-log relation varies from 
about 2 at the smallest masses to about 3 at 0.5-1 Mq. This allows us to calculate 
the critical descreening density (Eq. |20l) in terms of the parent star mass. The 
result is plotted for three different mass dependences of ^ {M,t) in Fig. [2] and for 
an assumed star- ISM velocity of 26 km s~^ 
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T3 



10^ I ' ' ' I 

^■^ stellar mass (MJ 

Figure 2: Critical interstellar density for descreening, defined as an event during 
which the astrosphere size is equal to the conventional liquid- water habitable-zone 
distance, as a function of parent-star mass in solar units. Three mass dependences 
of the stellar wind ram pressure are shown. Comparison with Fig. 1 shows that 
descreening densities are at or larger than the density above which gravitational 
focusing dominates. Even in the case in which low-mass stars have weaker winds 
((^ <x M), the descreening density still increases with decreasing mass, due to the 
severe dependence of the HZ distance on stellar mass. 
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A planetary system orbiting a star of mass 1 Mq will be descreened if it en- 
counters a region of density > 600 V20 cm^-^. On the other hand a system orbiting 
a star of mass 0.2 Mq must pass within a region of density > 10^ V20 cm^^ to be 
descreened, even if we assume no dependence of wind strength on mass or age. If 
we assume that the star- ISM velocity parameter, V20, scales with the kinematic ve- 
locity dispersion, then V20 could b e 2-3 times large r for bo th M stars (jpuchs et all 
200% and older solar-type stars (INordstrom et al.L |2004|) than for the Sun. Our 



results then imply that these stellar systems would require proportionally larger 
ISM densities to descreen the habitable zone. In this respect, our conclusions are 
somewhat conservative. 

Section [33] will argue that such dense regions will be rarely encountered. 

3.3 Frequency of Descreening as a Function of Stellar Mass 

For any reasonable characterization of the ISM, the frequency of encounters of a 
planetary system with a cloud of a given density decreases with increasing density. 
Unfortunately this statement is difficult to quantify because interstellar densities 
are normally very poorly determined. The densities are almost always derived 
from the observed column density of some tracer (e.g. dust, or a molecule like 
CO), a conversion factor from tracer to to total column density, and then a con- 
version to density by dividing the column density by a characteristic size of the 
region of interest. This latter step is especially uncertain because the apparent size 
is determined by the sensitivity limit of the observations (column density normally 
decreases outward on average) and by the ad opted di stance to the region. How- 



ever we can use the approach of Talbot and N ewman (|1977r) as applied to various 
catalogues of regions identified as discrete "clouds" in different surveys in order 
to see that, even given these uncertainties and selection effects, the frequency of 
occ urrence of dense ISM regio ns must de crease r apidly with increasing density. 



urrence or dense laM regio ns must de crease r apidly witn increasing density. 
Talbot and NewmanI (Il977h used the iLvndsl (Il962 h catalogue of about 1800 



"dark clouds,"Q for which Lynds provided "opacity classes" and angular areas 
Grouping these clouds into bins of a given column density range using the opac- 
ities, and using the mean size R{N) for each group, gives the number of clouds 
with density in each density range. Plotted as a histogram, this gives the num- 
ber of clouds per unit volume of space per unit density range, denoted nci{N), 
where N is the internal particle density. If 1{N) is the average number of clouds 



*Regions on photographic plates that appear darker than the surrounding area because of re- 
duction in the number of visible stars by dust extinction. 
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encountered of density {N,N + dN) per unit line of sight distance, then nci{N) = 
[nR{N)'^ 1{N)]^^. Since the frequency of encounters is proportional to 1{N)^^ 
this gives a frequency v <^ R{N)^ nci{N). 



We reexamined the cloud property table given by iTalbot and NewmanI (|1977|) . 
but we use only the data for the three highest opacity classes III Using the derived 
average densities in each class, we find the distribution nc\{N) N^^-^ and R(N) <x 
^-0.8^ giving a cumulative frequency of encounters v(> A^) <x N^^-^. This agrees 
with the exponent of —2.3 obtained using directly the /(> A^) values in their table, 
or the scaling of their quantity Ne which is the number of expected encounters by 
the Sun over 4.6 Gyr. 

This result is very uncertain because the Lynds cloud sample, besides its small 
size, suffers from a number of selection effects: 

1. Clouds with angular sizes less than about 2-5 arcmin cannot be detected 
on Palomar Schmidt plates because the number of background stars is too small. 
This corresponds to a minimum detectable cloud size that increases with distance. 
Since smaller clouds tend to be denser, the density probability distribution is prob- 
ably too flat because of this effect. 

2. Column density detection is somewhat restrictive, because regions with ex- 
tinction less than about one magnitude cannot be detected as a decrease in back- 
ground star density, while at extinctions above about 5-6 magnitudes clouds look 
black independent of their extinction. The latter effect causes the catalogue to 
lump clouds of densities above the saturation threshold into the bin of highest col- 
umn density, again making the apparent probability distribution of cloud densities 
appear too flat. 

3. The increasing number of foreground stars makes it progressively more dif- 
ficult to detect clouds farther than several hundred parsecs. This effect produces an 
underestimate of the extinction to distant large, and usually lower-density, clouds. 

4. Nearby clouds likely obscure more distant clouds, so th at some small clouds 



located in front of larger clouds are not detected, and so on (IScalo and Lazarianl . 



19961) . 

5. Distances to individual clouds detected in extinction are very uncertain. 
A way around this would be to examine only clouds in nearby complexes, like 
Taurus, which are all at about the same distance, but then the number of clouds 



would be small for statistical purposes (see IScalol |1985|) . iTalbot and Newman 



^The points for "Standard cloud," actually derived from reddening statistics, "CO" and "HI" 
are ill-defined averages for entire surveys, and the lower opacity class extinction clouds suffer 
from selection effects as explained below and as can be seen in Fig. 1 of Talbot and Newman. 
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1977h assumed they were all at the same distance. 



For these reasons we have also derived the relative frequency of encounters 
with regions above a given density using the FCRAO Outer Galaxy S urvey of car- 
bon m onoxide emission, which gives data for over 10,000 clouds (IHeyer et al. . 



20011). This survey is subject to selection effects analogous to 1 (angular reso- 
lution limit) and 2 (CO can only be excited above a density of about 100 cm^ 
but the [^^CO] line saturates at large column densities; also there is a sensitivity 
threshold on column density) above. The CO sample does not suffer from un- 
certainties 3-5 above because the clouds are detected using a spectral line, and 
the velocity information not only separates clouds along the same line of sight 
but allows distance estimates based on a kinematic model (rotation curve) for the 
Galaxy. We used the cloud properties from this survey (available electronically) 
to construct a density probability distribution nc\{N), but chose to include only 
clouds with distances less than 5 kpc, since a plot of the number of clouds vs. dis- 
tance indicated that incompleteness becomes serious beyond this point. We note 
that the clouds were identified by decomposition of larger structures, so the cata- 
logue does not contain larger, and on average less dense, regions within which the 
catalogued clouds were contained (Heyer, private communication). This should 
not affect our results much because these larger regions have average densities 
that are mostly below the critical descreening density for any but the most mas- 
sive parent stars considered here. Following Heyer et al., we take the calculated 
CO luminosity Leo to be proportional to the mass of the cloud and then divide 
by the size cubed in order to estimate the average internal density. Following this 
procedure, we find that nci{N) declines rapidly with density. If a power law is fit 
to the histogram, excluding the highest and lowest bins, the cumulative probabil- 
ity distribution function at the high-density end (2 x 10-^ to 2 x lO'^ cm^^) has a 
power-law slope of about —3, similar to the result found for the Lynds clouds. We 
also find the ^(A'^) °^ A^^^ for this sample; this result is in part a reflection of the 
same selection effect as for the Lynds cloud sample, since the surveys tend to be 
restricted in column density, yielding such a correlation. Nevertheless we adopt 
it as an average correlation. This gives a density dependence of the cumulative 
encounter rate of A'^^^, compared to A^^^-^ from Talbot and Newman. We also per- 
formed the same calculation but assuming all the clouds are at the same distance, 
to mimic the selection effect due to the unknown distances in Talbot and New- 
man's calculation, and find that nci(> A'^) becomes basically flat, but the distortion 
should be less severe for Talbot and Newman because the range in distances is 
smaller. 

We can also estimate the absolute frequency of descreening events for solar- 
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mass stars from the iHeyer et al\ (|2001h CO cloud catalogue and compare to Tal- 
bot and Newman's result for extinction clouds. For the 7900 CO clouds within 
5 kpc, we estimate a volume surveyed of 0.86 kpc-^, assuming the clouds are dis- 
tributed in a uniform disk with a half-thickness of 50 pc, which is close to most 
estimates for the scale height of the CO distribution perpendicular to the plane of 
the Galaxy. This yields the total cloud space density n^. Most of these clouds 
are at average internal densities above 300 cm^^ (the number declines rapidly 
below this density, suggesting selection effects), which is roughly the estimated 
descreening density for the solar system (Fig. [2l). The frequency of encounters is 
approximately v = 7r{R^)nc\v, where v is the mean star-cloud velocity (see Tal- 
bot and Newman, 1977)o We find a mean square cloud radius (R^) = 1.2 pc'^. 
Adopting a mean star-cloud velocity of 20 km s^^ gives a descreening frequency 
of 0.7 Gyr^. This is about six t imes smaller than ob t ained from including only 



the Lynds class 3-6 clouds from ITalbot and NewmanI (|1977L Table 1). Since we 



found that placing the Heyer et al. clouds all at the same distance significantly 
flattens the distribution, the assumption of a single distance could be part of 
the discrepancy. A major difference arises from the fact that the mean square sizes 
in the Lynds sample are much larger than found for the Heyer et al. sample. Such 
differences are perh aps expected beca use of the different selection criteria and se- 



lection effects (see IHeyer et a/.ll2001l for discussion in the CO sample), and the 
different distances to which the surveys extend, but this expectation gives us no 
basis to decide which is the more realistic estimate. The difference should prob- 
ably be considered a reflection of the inherent uncertainty involved in estimating 
the descreening frequency. Both results are lower limits because of the resolu- 
tion effect discussed earlier, by which smaller, denser clouds are missed at larger 
distances, but we do not know the size of the corresponding correction factor. 

Given these considerations, we think that the descreening frequency for the 
Sun cannot be estimated to better than an order of magnitude, but probably lies 
between 1 and 10 Gyr^^. Since the critical descreening density scales approxi- 



*Our frequency estimates assume a random distribution of clouds, but aren't much changed 
by the fact that these clouds, as well as those in Lynds (1962) sample, are actually clustered into 
larger structures. The mean encounter frequency is roughly the volume filling fraction of the 
clouds divided by their radius, and neither of these quantities is affected by clustering as long as 
the clouds don't overlap. The time series of encounters will be radically different in the clustered 
case, with long periods without descreening punctuated by clusters of frequent descreenings, but 
the average rate of encounters is similar if taken over a sufficiently long time interval (^ 1 Gyr 
here). For recent history, however, this clustering effect would be important. We discuss the 
details of the time history during intrusion into a large complex using hydrodynamic simulations 
in a separate publication. 
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mately as (Eq. [T6l) . and nc\{> N) decreases rapidly with A'^, we see that the 
descreening frequency declines extremely rapidly with decreasing stellar mass. 
Using the dark cloud nc\{> N) power-law slopes found above, we get a mass de- 
pendence of M^, while the CO cloud slopes yield M^. 

It should be pointed out that these results only apply to relatively dense clouds 
that can be seen in extinction or CO, and we do not know how the frequency scales 
with density for smaller densities. In addition, our estimate d rates for dense re- 
gions are probably underestimates, since both the Lynds and iHeyer et al\ (|200l|) 
cloud catalogues select against dense clouds because of resolution effects, as dis- 
cussed above. 

Since we showed earlier that an M star planet must encounter a region of 
density ~ 1 x 10^ cm^^ (^j. ^ O.SMq) to 3 x 10^ cm^^ (f^j. _ O.IM©) 
for descreening, while a solar-type star requires only ~ 600 cm^-^, the frequency 
estimates given above imply that descreening will occur 10^ to 10^ times less 
frequently for M star planets (i.e., never). Even for somewhat higher mass stars, 
the frequency of descreening is apparently such a strongly decreasing function 
of density that, given the crude estimate of 1-10 descreening encounters per Gyr 
for the Earth, even stars of spectral types late G and early K should be relatively 
immune from descreening. 



4 Summary and Conclusions 

We have used a ram-pressure balance model for calculating the size of astro- 
spheres in the apex direction around solar- and subsolar-mass stars to compute the 
minimum interstellar density required to "descreen" a planet in the habitable zone. 
We showed that gravitational focusing of incoming particle trajectories causes an 
enhancement of the ISM ram pressure and a consequent dependence of the as- 
trosphere size on the parent-star mass when the local ISM density rises above 
~ 100 M^^ cm^^, where M is the stellar mass in solar units. We have included 
the effects of gravitational focusing in calculating the ISM densities required to 
descreen planets in the habitable zone of solar-like and lower-mass stars and find 
a critical descreening density of roughly 600 M^^ cm^-^. Finally, we have esti- 
mated the frequency of descreening encounters as 1-10 Gyr^^ for solar-type stars 
and 10^-10^ times less frequent for M stars, so habitable planets around M stars 
will most likely never be exposed to the unfiltered ISM flow. 

We caution that these results are sensitive to the parameters of the stellar winds 
of low-mass stars. Studies of mass-loss rates for M stars have not yet well con- 
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strai ned the dependence of the mass-loss rate on mass or age. Preliminary results 
(e.g. IWood et alL I2002L l2005h suggest that the wind strength may increase with 
decreasing age and mass, but examples to the contrary do exist. While we neglect 
the age dependence of the wind strength for this discussion, we do find that even 
a mass dependence in which the mass-loss rate is linearly dependent on the stellar 
mass does not change our results (see Fig. 2). Deeper surveys of stellar mass- 
loss rates as a function of mass and age would be beneficial to understanding the 
astrospheric environment around low-mass stars, which are potenti ally the most 
common sites for habitable planets in the Galaxy dScalo et a/.Ll2007|) . 
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